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Abstract 

We calculate the density-dependent spin-asymmetry energy S(kf) of isospin-symmetric 
nuclear matter in the three-loop approximation of chiral perturbation theory. The interac- 
tion contributions to S(kj) originate from one-pion exchange, iterated one-pion exchange, 
and (irreducible) two-pion exchange with no, single, and double virtual A-isobar ex- 
citation. We find that the truncation to l-zr-exchange and iterated l-7r-exchange terms 
(which leads already to a good nuclear matter equation of state) is spin-unstable, since 
S(kfo) < 0. The inclusion of the chiral 7riVA-dynamics guarantees the spin-stability of 
nuclear matter. The corresponding spin-asymmetry energy S(kf) stays positive within a 
wide range of an undetermined short-range parameter S§ (which we also estimate from 
realistic NN-potentials). Our results reemphasize the important role played by two-pion 
exchange with virtual A-isobar excitation for the nuclear matter many-body problem. Its 
explicit inclusion is essential in order to obtain good bulk and single-particle properties. 



PACS: 12.38.Bx, 21.30.-x, 21.65.+f 



In recent years a novel approach to the nuclear matter problem based on effective field theory 
(in particular chiral perturbation theory) has emerged. The key element there is a separation 
of long- and short- distance dynamics and an ordering scheme in powers of small momenta. 
At nuclear matter saturation density po — 0.16 fm -3 the Fermi momentum kfo and the pion 
mass m n are comparable scales (fc/o — 2m,,-), and therefore pions must be included as explicit 
degrees of freedom in the description of the nuclear many-body dynamics. The contributions 
to the energy per particle E(kf) of isospin-symmetric (spin-saturated) nuclear matter as they 
originate from chiral pion-nucleon dynamics have been computed up to three-loop order in 
refs.fl, 2]. Both calculations are able to reproduce correctly the empirical saturation point of 
nuclear matter by adjusting one single parameter (either a coupling g + gi ~ 3.23 [1] or a cut- 
off A ~ 0.65 GeV [2]) related to unresolved short-distance dynamics. 1 The novel mechanism 
for saturation in these approaches is a repulsive contribution to the energy per particle E{kf) 
generated by Pauli-blocking in second order (iterated) one-pion exchange. As outlined in section 
2.5 of ref.[2] this mechanism becomes particularly transparent by taking the chiral limit m n = 0. 
In that case the interaction contributions to E{kf) are completely summarized by an attractive 
fcj-term and a repulsive /cy-term where the parameter-free prediction for the coefficient of the 
latter is very close to the one extracted from a realistic nuclear matter equation of state. 

In a recent work [3] we have extended the chiral approach to nuclear matter by including 
systematically the effects from two-pion exchange with single and double virtual A(1232)-isobar 
excitation. The physical motivation for such an extension is threefold. First, the spin-isospin- 
3/2 A(1232)-resonance is the most prominent feature of low-energy 7riV-scattering. Secondly, it 

1 Thc cut-off scale A serves the purpose to tunc the strength of an attractive zero-range NN-contact 
interaction. 



is well known that the two-pion exchange between nucleons with excitation of virtual A-isobars 
generates the needed isoscalar central NN-attraction [4] which in phenomenological one-boson 
exchange models is often simulated by a fictitious scalar "o""-meson exchange. Thirdly, the 
delta-nucleon mass splitting of A = 293 MeV is of the same size as the Fermi momentum 
kf ~ 2m w at nuclear matter saturation density and therefore pions and A-isobars should both 
be treated as explicit degrees of freedom. A large variety of nuclear matter properties has 
been investigated in this extended framework in ref.[3]. It has been found that the inclusion 
of the chiral 7riVA-dynamics is able to remove most of the shortcomings of previous chiral 
calculations of nuclear matter [2, 5, 6, 7]. For example, the momentum-dependence of the 
(real) single-particle potential U(p,kf) near the Fermi surface p — kf improves significantly. 
As a consequence of that the critical temperature of the first-order liquid-gas phase transition 
of isospin-symmetric nuclear matter gets lowered to the realistic value T c ~ 15 MeV. The 
isospin properties of nuclear matter improve also substantially when including the chiral 7riVA- 
dynamics. Instead of bending downward above po as in previous chiral calculations [2, 7], the 
energy per particle of pure neutron matter E n (k n ) and the (isospin) asymmetry energy A(kf) 
grow now monotonically with density (see Figs. 10,11 in ref. [3]). Good agreement with results of 
sophisticated many-body calculations and (semi)-empirical values has been found in the density 
regime p = 2p n < 0.2 fm~ 3 relevant for conventional nuclear physics. 

Given all that success of chiral perturbation theory for the nuclear matter many-body 
problem it is nevertheless still necessary to check the spin-stability of nuclear matter in that 
framework. Such an investigation is the subject of the present paper. We remind that in 
the past the criterion of spin-stability of infinite nuclear matter has required modifications 
of several phenomenologically very successful Skyrme forces [8, 9, 10]. For recent work on 
generalized symmetry energy coefficients in the context of phenomenological Skyrme forces, see 
also ref. [11]. 

Let us begin with defining the spin- asymmetry energy S(kf) of infinite nuclear matter. 
Consider this homogeneous many-nucleon system in a partially spin-polarized state such that 
the Fermi seas of the spin-up and spin-down nucleons are filled unequally high. With the help 
of the spin-projection operators (1 ± cr 3 )/2 such a spin-unsaturated configuration is realized by 
the substitution: 

0(k f -\p\) - l -^0{k,-\p\)+ l -^9{k l -\p\) ) (1) 

in the medium insertion. 2 Here, and k^ denote the different Fermi momenta of the spin-up 
and spin-down nucleons. Choosing = A;/(l ± t?) 1 / 3 (with 77 a small parameter) the total 
nucleon density p = (A; 3 + /c 3 )/3"7T 2 = 2kj/3n 2 stays constant. The expansion of the energy per 
particle of spin-polarized (isospin-symmetric) nuclear matter: 

E pol (k h fcj = E(kf) + V 2 S(k f ) + 0( V 4 ) , k u = k f (l± 77) 1 / 3 , (2) 

around the spin-saturation line (fcj = /q or 77 = 0) defines the spin-asymmetry energy S(kf). 
The obvious criterion for the spin-stability of nuclear matter is then the positivity of the spin- 
asymmetry energy, S(kf) > 0. The energy per particle at fixed nucleon density p must take on 
its absolute minimum value in the spin-saturated configuration. 

The first contribution to the spin-asymmetry energy S(kf) comes from the kinetic energy 
a/M 2 + p' 2 — M of a non-interacting relativistic Fermi gas: 

k 2 k A 

s{kf) = m'viM^ (3) 



2 Medium insertion is a technical notation for the difference between the in-medium and vacuum nucleon 
propagator. For further details, see section 2 in ref. [2]. 
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Figure 1: The two-loop one-pion exchange Fock diagram and the three- loop iterated one-pion 
exchange Hartree and Fock diagrams. The combinatoric factors of these diagrams are 1/2, 1/4 
and 1/4, in the order shown. Their isospin factors for isospin-symmetric nuclear matter are 6, 
12 and —6, respectively. 



with M = 939 MeV the (average) nucleon mass. The next term in this series, kj/16M 5 , is 
negligibly small at the densities of interest. 

Next, we come to interaction contributions to S(kf). The closed in-medium diagrams re- 
lated to one-pion exchange (Fock diagram) and iterated one-pion exchange (Hartree and Fock 
diagrams) are shown in Fig. 1. Differences in comparison to the calculation of the energy per 
particle E(kf) in ref. [2] occur only with respect to the factors emerging from the spin-traces over 
closed nucleon lines and the radii fcj^ = kf(l i??) 1 / 3 of the Fermi spheres to be integrated over. 
After some analytical calculation we find the following contribution to the spin-asymmetry 
energy S(kf) from the l7r-exchange Fock diagram in Fig. 1 (including its relativistic 1/M 2 - 
correction) : 



m 2 r5w 3 u 2 u , , . 9n 

arctan 2u + — 1 - 6m 2 In 1 + Au 2 ) 

.6 2 24 v 1 y ' 



M 2 



(4) 



Here, we have introduced the useful abbreviation u = kf/m n where m n = 135 MeV stands for 
the (neutral) pion mass. As usual f n = 92.4 MeV denotes the weak pion decay constant and we 
choose the value g& = 1-3 of the nucleon axial vector coupling constant in order to have a pion- 
nucleon coupling constant of g^N = 9aM/ f n = 13.2 [12]. In the second and third diagram in 
Fig. 1 the l7r-exchange interaction is iterated (once) with itself. These second order diagrams 
carry the large scale enhancement factor M (the nucleon mass). It stems from an energy 
denominator that is equal to a difference of small nucleon kinetic energies. With a medium 
insertion at each of two equally oriented nucleon propagators we obtain from the three-loop 
Hartree diagram in Fig. 1 the following contribution to the spin-asymmetry energy: 

S(kf) = ?r ^ Mm ; { ( I5u + —) ln(l + Au 2 ) - Uu - 16u 2 arctan 2^1 . (5) 
6(47r/ 7r ) 4 IV 2uJ J 

The right Fock diagram of iterated l7r-exchange (see Fig. 1) with two medium insertions on 
non-neighboring nucleon propagators gives rise on the other hand to a contribution to the 
spin-asymmetry energy of the form: 

... , 7rg\MmU32u 3 7u ( 11 2u\ , „ 2 . (\ rt \ . „ 2N 

r u arctan x — arctan 2x1 1 . . 

arctan 2u + I ax — > .(6) 

Jo u 1 + 2x 2 ) J J 



/32m 4 

+ f — 3 ) arctan u + (3 + 4m 2 ) 
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This expression does not include the contribution of a linear divergence / °° dl 1 of the momentum- 
space loop-integral. In dimensional regularization such a linear divergence is set to zero, whereas 
in cut-off regularization it is equal to a momentum space cut-off A. The additional term specific 
for cut-off regularization will be given in eq.(13). An in-medium diagram with three medium in- 
sertions represents Pauli-blocking effects in intermediate NN-states induced by the filled Fermi 
sea of nucleons. The unequal filling of the spin-up and spin-down Fermi seas shows its conse- 
quences in the spin-asymmetry energy. After some extensive algebraic manipulations we end 
up with the following double-integral representation of the contribution to the spin-asymmetry 
energy S(kf) from the Hartree diagram in Fig. 1 with three medium insertions: 



S(k f ) 



(4tt^)%3 



r-U r\ 

J dxx 2 J dy 



x 



x 



2s 2 + s 4 
1 + s 2 



-21n(l + s 2 ) 



2uxy(3u 2 - 5x 2 y 2 ) 

{v? — x 2 y 2 ) 

Au 2 H s 5 (4s' -3s) 

H + 

3(1 + s 2 ) 2 



{u 2 + 5x 2 y 2 )H 
2uxy + (u 2 — x 2 y 2 )H 

rA 



(5 + s 2 )(3s 2 - 8ss' + 8s' 2 ) + 8s(l + s 2 ){s" - 5s' + 3s) 



3(1 + s 2 ) 3 



where we have introduced several auxiliary functions: 

u + xy 



H = ln 



u — xy 



xy + \ju A — x A + x 2 y 2 , 



u 



as 

du 



u 



du 2 



(7) 



(8) 



Note that eq.(7) stems from a nine-dimensional principal- value integral over the product of 
three Fermi spheres of varying radii = kf(l ± rf) 1 ^ which has been differentiated twice 
with respect to rj at rj = 0. Of similar structure is the contribution to S(kf) from the iterated 
l7r-exchange Fock diagram with three medium insertions. Because of the two different pion 
propagators in the Fock diagram one ends up (partially) with a triple-integral representation 
for its contribution to the spin-asymmetry energy: 



S(k f ) 



g\Mm\ 
24(4vr^) 4 ^ 



£dx^G(3G 20 - 2G n + 3G 02 - 8G i - 3G) 



+3G 2 W -2G 01 G 10 + G 2 +Ax 2 



i dz yz6{y 2 + z 2 



l ^ J-i \yz\y/y 2 + z 2 — 1 



x 



x 



2sH 3 (8s't - 3st - As't') s 2 [t 2 - ln(l + t 2 )] 



+ 



(l + s 2 )(l + t 2 ) (1 + s 2 ) 2 

(3 + s 2 )(8ss' - 3s 2 - 8s /2 ) + 4s(l + s 2 )(6s' - 3s - 2s") 



(9) 



Here, we have split into "factorizable" and " non-factorizable" parts. These two pieces are 
distinguished by whether the (remaining) nucleon propagator in the three-loop Fock diagram 
can be canceled or not by terms from the product of 7rAMnteraction vertices. The factorizable 
terms can be expressed through the auxiliary function: 



G = u(l + u 2 + x 2 )-^[l + (u + x) 2 }[l + (u-x) 2 ] ln | + { M + ^ , 

Ax 1 + (u - x) 2 



(10) 



and its partial derivatives for which we have introduced a (short-hand) double-index notation 

, „■ d i+j G 



Gij = x l u 3 



dx l dui ' 



1 < i + j < 2 



(11) 
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Figure 2: The spin-asymmetry energy S(kf) of nuclear matter versus the nucleon density 
p = 2k^/3n 2 . The full line shows the result of a calculation up to fourth order in small 
momenta including l7r-exchange and iterated l7r-exchange. The cut-off scale A = 0.61 GeV has 
been adjusted to the saturation point: po — 0.173 fm -3 , E(kf ) = — 15.3 Me V. The isospin- 
asymmetry energy has the value A(kfo) = 38.9 MeV at saturation density. The negative values 
of S(kf) indicate the spin-instability of nuclear matter in this approximation. 



For the presentation of the nonfactorizable terms one needs also copies of the quantities s and 
s' defined in eq.(8) which depend (instead of y) on another directional cosine z: 

dt 

t = xz + Vu 2 - x 2 + x 2 z 2 , t' = u — . (12) 

au 

For the numerical evaluation of the dydz-dovhle integral in eq.(9) it is advantageous to first 
antisymmetrize the integrand both in y and z and then to substitute z = y/y' 2 C 2 + 1 ~ V 1 - 
This way the integration region becomes equal to the unit-square < y, £ < L In the chiral 
limit m n = the fourth order contributions in eqs.(5-9) sum up to a negative kj-term of 
the form: 5(^)1^=0 = (g A k f /4vr U) 4 (M/ '135) (16tt 2 + 273 - 936 In 2). The formulas for the 
contributions to the spin-asymmetry energy S(kf) written down here in eqs.(3-9) have a strong 
similarity with the analogous contributions to the isospin-asymmetry energy A(kf) collected 
in eqs. (20-26) of ref. [2]. This is not surprising in view of the spin-isospin structure a ■ Vr a of 
the pion-nucleon coupling. Finally, we give the expression for the linear divergence specific to 
cut-off regularization: 

2g\M_A l3 
(4ttA) 

to which only the iterated l7r-exchange Fock diagram (with two medium insertions) has con- 
tributed. In the case of the Hartree diagram the linear divergence drops out after taking the 
second derivative with respect to rj. 



S ( k f) = -7^ k h (13) 
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Now we can turn to numerical results. In Fig. 2 we show the spin-asymmetry energy S(kf) 
of nuclear matter as a function the nucleon density p = 2kj/3-ir 2 . The full line corresponds 
to a calculation up to fourth order in small momenta. It includes besides the kinetic energy 
term eq.(3) the contributions from static l7r-exchange and iterated l7r-exchange. For reasons of 
consistency we have dropped the small relativistic l/M 2 -correction in eq.(4) since it is of fifth 
order in the small momenta kf and m^. The cut-off scale A = 611 MeV has been adjusted to 
the nuclear matter saturation point p = 0.173 fm~ 3 and E(kf ) = —15.3 MeV [7]. The value of 
the nuclear matter compressibility, K = k^E'^kfo) = 252 MeV, as predicted in this framework 
is consistent with a recent extrapolation from giant monopole resonances of heavy nuclei in 
ref.[13], which gave K = (260 ± 10) MeV. One can read off from Fig. 2 a negative value of the 
spin-asymmetry energy at saturation density: S(kf ) = S'(2m 7r ) = —20.1 MeV. It reveals the 
spin-instability of nuclear matter in this restricted framework. The small positive values of S(kf) 
at very low densities p < 0.03 fm -3 result from the incomplete cancellation between the kinetic 
energy term eq.(3) proportional to k 2 and interaction contributions which grow in magnitude 
like /c 3 (or higher powers of kf). The largest negative contribution to S'(2m 7r ) = —20.1 MeV 
comes from the term linear in the cut-off A, eq.(13), which amounts to —35.5 MeV at saturation 
density kf = 2m 7T . Despite this overwhelming negative contribution to S(kf) one should not 
jump to the conclusion that the use of cut-off regularization were the reason for the spin- 
instability of nuclear matter in this calculation to fourth order. In essence the cut-off A is merely 
parameterizing the strength of a zero-range (S-wave) NN-contact interaction with a fixed ratio 
of its total-isospin 1 = and 1 = 1 components. The isospin factor (ti • r 2 ) 2 = 3 — 2t\ ■ r 2 of the 
iterated l7r-exchange diagram determines this ratio as 9 : 1. The quality of such a fixed ratio 
can be tested via the isospin-asymmetry energy A(kf) (for explicit expressions see section 3 in 
ref.[2]). The predicted value A(kfo) = 38.9 MeV for the isospin-asymmetry energy at saturation 
density of the present fourth order calculation employing cut-off regularization agrees within 
10% with a recent empirical determination, A(kf ) = (34 ± 2) MeV [13]. The prescribed ratio 
9 : 1 between the 1 = and 1 = 1 components of the emerging contact-interaction works 
therefore rather well. A key observation in this context is that terms linear in the density p 
satisfy the relation: 3S(kf)\[ n + 3A(kf)n n + 2E(kf)n n = 0. From that one can conclude that 
even with two free parameters in the contact-interaction (adjusted to the empirical values of 
E(kfo) and A(kf )) nuclear matter will be spin-unstable in a calculation restricted to fourth 
order in small momenta. For example, in dimensional regularization where the linear divergence 
Jq 00 dl 1 is set to zero one would instead employ a two-parameter momentum-independent NN- 
contact interaction. The adjustment to the empirical nuclear matter saturation point E(kfo = 
2m n ) = —15.3 MeV introduces then the explicit term linear in density £'(2m 7r )ii n = —177 '.4 MeV 
while the constraint from the empirical isospin-asymmetry energy A(kf ) = 34 MeV requires 
in addition A(2m 7T )\ in = 148.9 MeV. These fixed contributions imply via the abovementioned 
relation a term linear in density for the spin-asymmetry energy of S , (2m 7r )i in = —30.6 MeV. It 
is almost as negative as the value —35.5 MeV obtained in cutoff regularization with only one 
adjusted cut-off scale A and therefore the spin- instability of nuclear matter would still remain. 

The foregoing discussion has shown that it is mandatory to include contributions of fifth 
order in the small momentum expansion in order to be able to reach the spin-stability of 
nuclear matter. In the three-loop approximation the terms of fifth order are generated by 
(irreducible) two-pion-exchange between nucleons. The corresponding one-loop diagrams for 
elastic NN-scattering are shown in Fig. 3. Since we are counting the delta- nucleon mass splitting 
A = 293 MeV (together with kf and m n ) as a small momentum scale the diagrams with single 
and double virtual A(1232)-isobar excitation shown in Fig. 4 belong to the same order. The 
non-relativistic nucleon and delta propagators are both counted in this scheme as an inverse 
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Figure 3: One- loop diagrams of irreducible two-pion exchange between nucleons. 
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Figure 4: One-loop two-pion exchange diagrams with single and double A(1232)-isobar excita- 
tion. Diagrams for which the role of both nucleons is interchanged are not shown. 



small momentum. 

By closing the two open nucleon lines of the one- loop diagrams in Figs. 3,4 to either two 
or one ring one gets (in diagrammatic representation) the Hartree or Fock contribution to the 
energy density of nuclear matter. The Hartree contribution to the spin-asymmetry energy 
S(kf) vanishes identically because the relevant 27r-exchange NN T-matrix in forward direction 
is spin-independent [4, 14]. The Fock contribution on the other hand is obtained by integrating 
the spin- and isospin-contracted T-matrix over the product of two Fermi spheres of radii k^i = 
kf(l i??) 1 / 3 . We separate regularization dependent short-range contributions to the T-matrix 
(originating from the ultra-violet divergences of the one-loop diagrams in Figs. 3,4) from the 
unique long-range terms with the help of a twice-subtracted dispersion relation. The occurring 
subtraction constants give rise to a contribution to the spin-asymmetry energy of the form: 



S(k f ) = (B 3 - 6B 



n,3 



i.3 u5 
K f i C _f_ 

3M 2 5 M 4 



(14) 



The two dimensionless parameters B% = —7.99 and B n ^ = —0.95 have been adjusted in 
ref.[3] to few empirical nuclear matter properties (such as the maximal binding energy per 
particle —E(kf ) = 16MeV and the isospin-asymmetry energy A(kf ) = 34MeV with kf = 
261.6 MeV). Again, we recognize in the first part of eq.(14) the relation 3S(kf)\ in + 3A(kf)\- m + 
2E(kf)\i n = for terms linear in density p = 2k^/3ir 2 . The other subtraction constant S5 in 
front of the kj/M 4 -teim is (a priori) not constrained by any property of spin-saturated nuclear 
matter. The long-range parts of the 27r-exchange (two-body) Fock diagrams can be expressed 
as a dispersion-integral: 



Qk 3 f lQk 5 f 3u 3 , / Akl^ 



+lm(V T + 3W T ) 



6^ 3 k f + 4fik 3 f 



2k 



/// 



(15) 
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Figure 5: Hartree and Fock three-body diagrams related to 27r-exchange with single virtual 
A-isobar excitation. They represent interactions between three nucleons in the Fermi sea. In 
the case of isosp in-symmetric nuclear matter the isospin factors of these diagrams are 8, 0, and 
8, in the order shown. The combinatoric factor is 1 for each diagram. 



where ImVc, ImW / c', ImVr and ImWr are the spectral functions of the isoscalar and isovector 
central and tensor NN-amplitudes, respectively. Explicit expressions of these imaginary parts 
for the contributions of the triangle diagram with single A-excitation and the box diagrams 
with single and double A-excitation can be easily constructed from the analytical formulas given 
in section 3 of ref.[4]. The //- and /^-dependent weighting functions in eq.(15) take care that 
at low and moderate densities this spectral-integral is dominated by low invariant 7T7r-masses 
2m n < fi < 1 GeV. The contributions to the spin-asymmetry energy S(kf) from irreducible 
27r-exchange (with only nucleon intermediate states, see Fig. 3) can also be cast into the form 
eq.(15). The corresponding non- vanishing spectral functions read [14]: 



lmWc(ifj) 



Am? 



3vr/i(4A) 



4ml(l+4gl-5g\)+f,\23g\- 
6g\J fi 2 — 4ml 

ImVr(?/i) = - 



10^ 



1) + 



48g\ml 



ji 2 — 4m\ 



(16) 



(17) 



Next, we come to the additional 27r-exchange three-body terms which arise from Pauli 
blocking of intermediate nucleon states (i.e. from the —2ird(k^^ — \ f>\) terms in the in-medium 
nucleon propagators [2]). The corresponding closed Hartree and Fock diagrams with single 
virtual A-excitation are shown in Fig. 5. For isospin-symmetric nuclear matter their isospin 
factors are 8, 0, and 8, in the order shown. The contribution of the left three-body Hartree 
diagram to the spin-asymmetry energy S(kf) has the following analytical form: 



S(k f ) 



Q \ ii 2 

- + 4u 2 j ln(l + 4u 2 ) - 2u\l + 3() - 8m 2 - -^-^ 



(18) 



The delta propagator shows up in this expression merely via the (reciprocal) mass-splitting 
A = 293 MeV. Furthermore, we have already inserted in eq.(18) the empirically well-satisfied 
relation g w NA = ^9-kn/V2 for the 7riVA-coupling constant. The parameter ( = —3/4 has 
been introduced in section 2 of ref.[3] in order to reduce a too strongly repulsive p 2 -term 
in the energy particle E(kf). It controls the strength of a three-nucleon contact interaction 
~ ((g\/Af%) (NN) 3 which has the interesting property that it contributes equally but with 
opposite sign to the energy per particle E(kf) and the spin-asymmetry energy S(kf). The 
contribution of the right three-body Fock diagram in Fig. 5 to the spin-asymmetry energy 
S(kf) can be represented as: 



S(k f ) 



4 6 
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Figure 6: The spin-asymmetry energy S(kf) of nuclear matter versus the nucleon density 
p = 2k 3 f /3-K 2 . In comparison to Fig. 2 the effects from 27r-exchange with single and double virtual 
A-isobar excitation are now included. The full, dashed and dashed-dotted curves correspond 
to the choices S$ = 0, —7.45 and —30 of the short-range parameter S§. The positive values of 
S(kf) ensure the spin-stability of nuclear matter. 



+2Gs{3G$ + 8G501 — 3G502 + 2G511 — 3G520) + 2GtoiGtio 

— 7G^ 01 — 3Gy 10 + Gt(3Gt + 8Gtoi — 3GVo2 + 2Cni — 3Gt2o)} , 

with the two auxiliary functions: 

4ux 



G s = 



3 

+ (x 2 



(2u — 3) + Ax arctan(w + x) + arctan(« 



x 



u 



l)ln 



1 + (u + xf 
1 + (u — x) 2 ' 



(19) 



(20) 



Gt — 



— (8u 2 + 3x 2 )-— (l + u 2 ) 2 



2\3 



U ) 



X/ 



2x 

x 4 + (l-3u 2 )(l+u 2 -x 2 ) 



In 



l + (u + x) 
1 + (u — x) 



(21) 



The double-indices on Gj and Gt have the same meaning as explained in eq.(ll) for the 
function G. It is remarkable that the three-body Hartree and Fock diagrams with virtual A- 
excitation (shown in Fig. 5) lead both to identical contributions to the spin-asymmetry energy 
S(kf) and the isospin-asymmetry energy A{kf) (compare eqs.(18,19) here with eqs. (48,49) in 
ref.[3]). This feature originates from the specific spin-isospin structure of the (non-relativistic) 
7r./VA-transition vertex [4]. 

In Fig. 6 we show again the spin-asymmetry energy S(kf) of nuclear matter as a function of 
the nucleon density p = 2kl/3ii 2 . The full line includes all the contributions from chiral In- and 
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27r-exchange written down in eqs. (3-9, 14-19). The (yet undetermined) short-range parameter 
S 5 has been set to zero, S 5 = 0. We note as an aside that the term linear in the density and the 
cut-off A eq.(13) is of course now not counted extra since the parameter P 3 — 6B n 3 = —2.29 
[3] collects all such possible terms. One observes in Fig. 6 a positive spin-asymmetry energy 
S(kf) which rises monotonically with the density p. The inclusion of the chiral 7riVA-dynamics 
has made nuclear matter spin-stable. This result is of utmost importance since it assures that 
the improved nuclear matter properties reported in ref. [3] are based on a spin-stable ground 
state. It is also interesting to look at numerical values of S(kf) and their decomposition. At a 
Fermi momentum of kf = 2m^ (corresponding to p = 0.173 fm -3 ) the spin-asymmetry energy 
is now S(2m n ) = 67.3 MeV (setting S 5 = 0). The most significant changes in comparison to 
the previous fourth order calculation come from the term linear in density which gets reduced 
(by a factor of 2.1) to —17.0 MeV and the two-body Fock and three-body Hartree contributions 
eqs.(15,18) which amount together to 54.6 MeV + 13.7 MeV = 68.3 MeV. We note also that 
about one quarter thereof (namely 16.2 MeV) stem from the three-body contact interaction 
proportional to ( = — 3/4. 

The remaining open question concerns the size of the short-distance parameter S5 in eq.(14). 
Large negative values could again endanger the spin-stability of nuclear matter. In order to get 
an estimate of 65 we bring into play the complete set of four-nucleon contact-couplings written 
down in eqs. (3,4) of ref. [15]. This set represents the most general short-range NN-interaction 
quadratic in momenta and it involves seven low-energy constants C±, . . . , C7. After computing 
the spin-asymmetry energy S(kf) from the corresponding contact-potential in Hartree- Fock 
approximation we find: 

S 5 = ^(C 2 -4C 1 + 3C 4 + C 7 ) 

M 4 r 1 

= — [2C(^i) - 2C( 1 S ) - C( 1 Pi) + C( 3 P ) + 3C( 3 Pi) + 5C( 3 P 2 )] . (22) 

In the second line of eq.(22) we have reexpressed the relevant linear combination of 612,4,7 
through the so-called spectroscopic low-energy constants which characterize the short-range 
part of the NN-potential in the spin-singlet and spin-triplet S- and P-wave states. In that 
representation we obtain from the entries of table IV in ref. [15] for the three high-precision 
NN-potentials CD-Bonn [16], Nijm-II [17] and AV-18 [18] the numbers S 5 = -7.10, -7.98 and 
—7.27. The dashed line in Fig. 6 shows the spin-asymmetry energy S(kf) which results from 
taking their average value S 5 = —7.45. The effective reduction of the spin-asymmetry energy 
by about 20% is harmless and the spin-stability of nuclear matter is equally well guaranteed. 
The dashed-dotted curve in Fig. 6 corresponds to the extreme choice S5 = —30. One can see 
that even with such a large negative S^-value, which exceeds the estimate from realistic NN- 
potentials by a factor of 4, the spin-stability of nuclear matter still remains preserved. We can 
therefore conclude that with the inclusion of the chiral 7rA^A-dynamics spin-stability becomes 
a robust property of nuclear matter. This is an important finding and it goes conform with the 
improved isospin properties observed in ref. [3]. 

In summary we have investigated in this work the spin-stability of nuclear matter in the 
framework of chiral perturbation theory. For that purpose we have calculated the density- 
dependent spin-asymmetry energy S(kf) of isospin-symmetric nuclear matter to three-loop 
order. The interaction contributions to S(kf) originate from one-pion exchange, iterated one- 
pion exchange, and (irreducible) two-pion exchange with no, single, and double virtual A- 
isobar excitation. We have found that the truncation to l7r-exchange and iterated l7r-exchange 
terms (which leads already to a good nuclear matter equation of state) is spin-unstable, since 
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S(kf ) < 0. This statement holds independently of the regularization scheme if the contact 
terms (generating contributions linear in density) are consistent with the empirical nuclear 
matter bulk properties: E(kfo) ~ — 16MeV and A(kf ) ~ 34MeV. The inclusion of the chi- 
ral 7rVA-dynamics on the other hand guarantees the spin-stability of nuclear matter. The 
corresponding spin-asymmetry energy S(kf) stays positive within a wide range of an undeter- 
mined short-range parameter S5, which we have also estimated from realistic NN-potentials. 
Our results reemphasize the important role played by two-pion exchange with virtual A-isobar 
excitation for the nuclear matter many-body problem. 
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